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We present two different strategies for developing a quantum information science platform, based 
on our experimental results with magnetic microtrap arrays on a magnetic-film atom chip. The 
first strategy aims for mesoscopic ensemble qubits in a lattice of ^ 5 /xm period, so that qubits can 
be individually addressed and interactions can be mediated by Rydberg excitations. The second 
strategy aims for direct quantum simulators using sub-optical lattices of ~ 100 nm period. These 
would allow the realization of condensed matter inspired quantum many-body systems, such as 
Hubbard models in new parameter regimes. The two approaches raise quite different issues, some 
of which are identified and discussed. 



I. INTRODUCTION 

Atom chips are a promising technology for develop- 
ing a novel QI science platform pQ. They combine the 
best of two worlds, on the one hand neutral atoms with 
their weak coupling to the environment and concomitant 
long coherence times [2[ |3] , and on the other hand the 
compactness and large-scale integration possibilities of 
solid-state like systems. 

Here we present our experiments using an array of in- 
dividually addressable magnetic microtraps defined on a 
magnetic-film atom chip |4j [5] . We discuss the scaling 
issues, and prospects to develop a quantum simulator, 
which for the first time is becoming technologically re- 
alizable. A series of theoretical proposals [6 - 8 has sug- 
gested that experimental techniques available today can 
be synthesized into an universal device capable of mim- 
icking the behavior of complex, many-spin systems. 

Typically, simulating quantum many-body physics has 
required a lattice structure into which individual parti- 
cles, whether neutral atoms, ions, or electrons, can be 
loaded and then coherently controlled. Analog quantum 
simulators can be used to find the ground state of in- 
teracting quantum many-body systems such as Hubbard 
models. Coherent control can also be used to prepare 
the simulator in an analogous initial state, after which 
its subsequent dynamic evolution naturally mimics that 
of the system of interest. Of greater versatility is the 
digital quantum simulator, in which quantum gates are 
operated to produce unitary time evolutions. This is also 
known as an universal quantum computer [9|. 

Both types of simulators can be implemented with 
neutral atoms in a lattice. For atom chip devices, lat- 
tices formed by magnetic potentials offer some promising 
technical prospects compared to other approaches, such 
as optical lattices. Arbitrary lattice configurations and 
spacings can be implemented without the need to con- 
sider the wavelength of the light or optical access close to 
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the chip surface. This becomes particularly critical when 
considering the exciting regime of sub-optical wavelength 
lattice spacings on the order of 100 nm. 

As we have recently demonstrated a working shift reg- 
ister on a magnetic lattice of 20 /mi spacing [4], the major 
questions we are addressing in this article can be summa- 
rized as "the physics of scaling down". We present our 
estimates for the required specifications of a magnetic- 
lattice-based register at 5 fim and, ultimately, 100 nm. 
We will be considering issues in the context of both the 
implementation of a digital quantum simulator composed 
of Rydberg gates and a true analog quantum simulator 
at a lattice spacing of 100 nm. 

For the micron-scale array of mesoscopic ensemble 
qubits, we describe our results on three-body decay lead- 
ing to sub-Poisson atom number fluctuations [10 . We 
also discuss our progress on sensitive detection of small 
atom numbers and the investigation of Rydberg exci- 
tation close to the chip surface. First experiments on 
Rydberg atoms on chips have already been performed, 
in which we measured surface-based electric fields using 
electromagnetically-induced transparency [TT] . 

The sub-optical scale lattices can potentially yield a 
novel realization of a Hubbard model system that can 
open up new parameter regimes. We discuss the scaling 
of the Hubbard model parameters with our atom chip in 
mind. We calculate the influence of the Van der Waals 
surface attraction and show that the magnetic trapping 
potential is strong enough to overcome this. We also 
calculate rates for tunneling into the surface, and show 
that for practical cases tunneling into the surface can 
be neglected. Finally, we estimate trap loss rates due 
to Johnson-noise-driven spin flips and explore ways to 
mitigate this effect. 

A non-trivial technical challenge which accompanies 
increasing minaturization is the question of effective de- 
tection of closely spaced atoms on the single-particle 
level. Although much of our research efforts have been to- 
wards engineering magnetic trapping potentials, we shall 
briefly discuss an approach to tailor electric fields, based 
on ideas inspired by earlier work on field emission from 
arrays of silicon tips [12] and from individual carbon nan- 
otubes [13]. This offers several opportunities, in partic- 
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ular for detection based on site-selective field ionization 
on a magnetic atom chip. 

II. MICROTRAP ARRAY ON A 
MAGNETIC-FILM ATOM CHIP 




FIG. 1: Upper: scanning electron microscopy image of the 
patterned magnetic film. The light grey areas correspond to 
the FePt pattern and the darker regions are the Si substrate. 
The size of the horizontal edges is 20 fim. Lower: schematic 
representation of the microtrap array on the magnetic-film 
atom chip. Small ultracold (~ /xK) clouds of Rb atoms are 
trapped above edges of the magnetized FePt-film pattern, at 
c ~ 10 fim from the surface. Lattice spacings are a = 22 /xm 
and b = 36 /im. 

We briefly review our previously reported results on an 
array of microscopic Ioffe-Prit chard type magnetic traps, 
defined by patterning a permanently magnetized film on 
a chip [4j [5] . The film with a thickness of 300 nm consists 
of FePt with a remanent magnetization of 670 kA/m ori- 
ented perpendicular to its surface, and is deposited on a 
300 /im-thick silicon substrate. Such high remanent mag- 
netization is an important advantage of FePt, in addition 
to its high coercivity, a high Curie temperature that al- 
lows mild baking of the vacuum system, and its corro- 
sion resistance. FePt is also very suitable for patterning 
via conventional lithographic techniques. Transfer of the 
pattern was made by UV lithography and plasma etching, 
resulting in the pattern shown in Fig. [I] after which the 
chip is coated with a reflecting gold layer. This process 
yields a two-dimensional array of 1250 traps/mm 2 . Each 
trap hovers 10-15 /im above the chip surface, see also 
Fig. [I] the traps are tightly confining, with trapping fre- 
quencies of ~ 10 kHz for 87 Rb atoms in the \F = mp = 2) 
state. 

The chip is mounted on a copper surface with an in- 
tegrated Z-wire, which is used for loading in conjunction 



with an external bias field. Radio frequency evaporative 
cooling is used to lower the temperature of the atoms to 
quantum degeneracy in the micro-Kelvin regime. Fig- 
ure]^ shows an averaged absorption image of the loaded 
magnetic lattice. Here, ~ 500 traps are loaded with be- 
tween 200 and 2500 atoms each. Individual traps are re- 
solved, with an optical resolution of ~ 7.5 jam (Rayleigh 
criterion), enabling individual detection of the number 
of atoms in each ensemble. Our microtrap array of near- 
degenerate mesoscopic atomic ensembles bridges the gap 
between optical lattices with a single or few atoms per 
site on the one end, and single macroscopic Bose-Einstein 
condensates with 10 3 — 10 8 atoms on the other end. 




200 nm 



FIG. 2: Absorption image of the loaded lattice, showing ~ 500 
traps loaded with 200-2500 atoms each [4]. 

It is possible to address individual sites and selectively 
empty one site using a focused laser pulse. The optical 
pumping transition F = 2 — >> F' = 2 of the D2 line is 
used to flip the spin of the atoms in a selected site into a 
magnetically anti-trapped state. This can be combined 
with parallel transport of atoms along the array, which 
has also been demonstrated [4] . External magnetic fields 
were dynamically varied to shift the potential minima 
of the micro-traps across the surface of the chip. Fig- 
ure [3] shows shifting over two lattice periods, using a va- 
cancy made by the addressing laser as a marker. During 
the shifting, the empty sites do not refill, demonstrat- 
ing there is no interaction between neighbouring sites. 
Atoms at « 10 fiK were transported 360 fim using this 
method, without significant additional heating. 




FIG. 3: Optical addressing of single lattice sites and atomic 
shift register, (a) Close up absorption image of the central 
region of the lattice. The trap at the focus point of the ad- 
dressing laser is circled in red. (b) A single site has been 
emptied by a 1 ms laser pulse, (c-d) The entire lattice is 
shifted two cycles by rotating the external magnetic bias field 

®- 

The magnetic microtrap array is an excellent starting 
point to develop a scalable quantum information science 
platform. In the remainder of the paper we discuss some 
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aspects of two different strategies that we are pursuing 
to achieve this goal. 

The first involves the redesign and scaling down of the 
array to arrive at square and/or triangular lattices of 
~ 5 jam period. This would allow the definition of an en- 
semble qubit at each lattice site, for which the transition 
between the logical states |0), |1) is a collective, single- 
atom excitation to another ground hyperfine state. This 
approach is based on proposals relying on dipole blockade 
between highly excited Rydberg atoms [6HS1 HU US] • 

The second strategy involves scaling down the array 
to sub-optical length scales of ~ 100 nm period. This 
may open up new parameter regimes in Hubbard model 
physics that are as yet out of reach for optical lattices, 
allowing the realization of direct quantum simulators of, 
for example, condensed-matter-inspired systems. In this 
approach the atoms will be trapped much closer to the 
surface. This makes the use of Rydberg atoms unrealistic 
and makes it necessary to carefully consider the atom- 
surface interaction. 



III. MESOSCOPIC ENSEMBLES FOR 
QUANTUM INFORMATION 

The use of mesoscopic ensembles for quantum infor- 
mation hinges on the use of strong, long-range interac- 
tions among highly excited Rydberg atoms to orches- 
trate switchable gates [6] [TJ [14]. In the dipole block- 
ade effect, the strong electric dipole moment or polar- 
izability of an atom excited to a Rydberg state shifts 
the energy levels of nearby Rydberg atoms, prohibiting 
the excitation of more than one atom to the Rydberg 
state within one so-called blockade radius. The strong, 
long-range, and switchable interaction fuels the interest 
in Rydberg atoms for implementing quantum gate pro- 
tocols and multi-particle entanglement [SJ [I5j [16] . 

Dipole blockade has been proposed as a way to define 
ensemble qubits and to create entangled mesoscopic en- 
sembles. So far the dipole blockade has been observed 
either in clouds much larger than a blockade radius or 
among separate single atoms [T71H9] . Recent experi- 
ments have demonstrated dipole blockade between sin- 
gle atoms up to 10 jam apart and CNOT gate opera- 
tion between two individually addressed neutral atoms. 
Resonant energy transfer has been demonstrated over 
40 jam [20]. For comparison, an atom in the ground vi- 
brational level of a 10 kHz microtrap is confined to less 
than 100 nm. Even an ensemble at a temperature of 
T = 2 jaK would have a rms size of only 200 nm. Thus 
the separation between mesoscopic ensembles in strongly 
confining microtraps are well within the Rydberg inter- 
action radius, making them a promising candidate qubit. 

In ensembles, qubits could be encoded in collective ex- 
citations produced using intratrap dipole blockade. We 
consider two trappable hyperfine ground states, for ex- 
ample the \F,m F ) = 1 1,-1) \A) and |2,1) \B) states 
in 87 Rb, which have nearly the same magnetic moment. 



At the special magnetic field value of 3.23 G these two 
states serve as "clock states" with vanishing first order 
differential Zeeman shift [2 . The logical qubit states 
in an ensemble of N atoms could be defined as follows: 
|0) = \A)® N and |1) = \A)^ N ~^\B). Thus the logical 
|1) is a single, collective hyperfine ground state excita- 
tion, which can be accessed via an intermediate collec- 
tive Rydberg excitation if the entire ensemble is within 
one blockade radius. The strong confinement in magnetic 
microtraps brings us in this interesting regime where the 
dipole blockade allows the creation of a single collective 
Rydberg excitation \A)®( N ~ 1 " ) |r), which can be mapped 
back on |1) by a 7T pulse [21] E2]. 

The prospect of using controlled Rydberg interactions 
also among neighboring microtraps have led us to re- 
design the array. Where the previous design had a par- 
allelogram unit cell with lattice spacings of 22 jam and 
36 /im, respectively (Fig. [T]), the new design features a 
square and a triangular lattice both with a period of 5 jam 
in all directions. Equal lattice spacings are a require- 
ment for truly 2D nearest-neighbour interactions. The 
5 jam spacing is chosen to balance strong Rydberg dipole- 
dipole interactions with sufficient optical resolution to re- 
solve individual traps with in- vacuum optics. Thus the 
high trap frequencies (> 10 kHz) and relatively small 
lattice constants of approximately 5 jam will allow us to 
investigate both intra- and inter-trap interactions. The 
new magnetic lattice patterns have been designed using 
a numerical optimization algorithm [23]. The resulting 
magnetization pattern for the square lattice is shown in 
Fig. [4^ a) in a downscaled version. 

The approaches sketched here raise several issues that 
we will discuss in more detail. 



A. Sub-Poisson atom number fluctuations 

In collective excitation schemes, an ensemble contain- 
ing TV atoms would benefit from a y/N enhancement 
in Rabi frequency, increasing the number of operations 
achievable within the coherence time [7 . However, the 
TV-dependence of the Rabi frequency adversely affects the 
fidelity of operations if the number of atoms fluctuates 
from shot to shot and/or from site to site. Intrinsic atom 
number fluctuations would typically lead to a Poisson 
distribution at best; however, the fluctuations can be re- 
duced to sub-Poissonian levels by three-body losses. In 
our case we find that three-body loss is the dominant trap 
loss mechanism due to the high trap frequencies and cor- 
responding high densities [ID] . 

We find theoretically and experimentally that the nor- 
malized atom number variance, or Fano factor, F = 
(N 2 )/(N) is reduced to the value 0.53(22), significantly 
below the value of 1 corresponding to Poisson statistics. 
Furthermore, the memory of any initial fluctuations is 
quicky erased. The Fano factor decays as the fifth power 
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of the remaining fractional atom number 77, 

F( V ) = l+ V 5 (Fo-iy (1) 

This means in practice that once half the atoms are lost, 
the Fano factor is usually below 1 and already close to 
its asymptotic value of 3/5. 

The occurrence of sub-Poissonian atom number fluc- 
tuations can be seen as an atomic analog to intensity 
squeezing in optics. By trapping a large number of dense 
mesoscopic ensembles in a lattice of microtraps which un- 
dergo rapid three-body decay, it was shown that three- 
body loss can be used to prepare small and well-defined 
numbers of atoms in each trap. Through sensitive ab- 
sorption imaging the shot-to-shot distribution of atom 
numbers was measured and found to have sub-Poissonian 
statistics for between 50 and 300 atoms per trap, in good 
agreement with a model for stochastic three-body loss 
which takes into account the fluctuations. 



B. Rydberg-surface interaction 

Little has so far been known about the influence of 
the nearby surface on Rydberg atoms created on an 
atom chip [24] [25] . We investigated such surface ef- 
fects [11 using excited-state electromagnetically induced 
transparency (EIT) [26, 27 and extending this technique 
to obtain EIT spectra at distances between « 20 — 200/mi 
from the chip surface. The position and width of the nar- 
row transmission resonance reflect the energy and life- 
time of the Rydberg state under investigation and pro- 
vide a sensitive probe of the atom-surface interaction. 

We find strong shifts of the Rydberg states of several 
10s of MHz for principal quantum numbers n « 30 at 
these distances. The magnitude and sign of the shift are 
determined by the polarizability of the Rydberg state; 
they are caused by a dipolar electric field created by a 
patch of Rb adsorbed on the chip which has been de- 
posited there over many experimental cycles during reg- 
ular operation of the experiment [28]. We also find that 
there is no broadening of the Rydberg resonances above 
the level given by the linewidth of the lasers used in the 
experiment, either by this patch of adatoms or by other 
sources. 

The surface effects due to adatoms could be prevented 
in future experiments by incorporating a magnetic field 
gradient to push the atoms away from the surface at the 
end of each experimental cycle, or the adatoms could be 
removed using e.g. light induced desorption [29]. Fur- 
thermore, use of other coating materials on the chip sur- 
face could decrease the dipole moment of the adatoms 
or increase the desorption rate. Finally, one could adjust 
the frequency of the coupling laser to compensate for any 
residual effects. 

The absence of any broadening and the proposed mea- 
sures for dealing with the level shifts caused by adatoms 
lead to the conclusion that surface fields do not cause any 



serious obstacles to creating Rydberg atoms at distances 
of order 10 /am from the surface of an atom chip. 



C. Sensitive detection of small atom numbers 

An important technical challenge lies in the state- 
resolved detection of either a single atom or a single 
collective hyperfine excitation of an ensemble. The di- 
rect detection of atoms in their Rydberg state would re- 
quire additional technology such as ion detection, etc. 
Although the Rydberg atom can be mapped back on a 
detectable ground state such as the \F = 2) by a 7r-pulse 
in the regime of active dipole blockade, this would still 
require the capability to detect at the single-atom level. 

Our group has recently demonstrated sensitive detec- 
tion of small (~ 10) numbers of atoms by reflective 
absorption imaging, using advanced image processing 
[30] . First, a fringe removal algorithm was used to re- 
duce imaging noise to the fundamental photon-shot-noise 
level. A maximum-likelihood estimator (MLE) was then 
used for optimal atom-number estimation. The MLE is 
optimal in the sense that it achieves the Cramer-Rao 
bound: a lower bound for the variance of any parame- 
ter estimate, independent of the exact procedure used to 
extract the information. In particular the variance in the 
estimated atom number was shown to be lower than one 
would obtain from a simple integration over a region of 
interest in the image. Both the detection sensitivity and 
the spatial resolution that were achieved were limited by 
the low (0.1) numerical aperture (NA) of the imaging op- 
tics. Increasing the NA is straightforward by mounting 
a high NA lens inside the vacuum chamber. Our esti- 
mate shows that single-atom, single-shot readout sensi- 
tivity using absorption imaging is entirely realistic. 

This capability will be important, apart from the state- 
selective readout of ensemble qubits, to verify the phe- 
nomenon of dipole blockade. Two important milestones 
will be to detect the blockade by one Rydberg atom of 
further excitation to the Rydberg level of all other atoms 
in the same potential well, followed by the detection of 
blockade of Rydberg excitation in adjacent wells. 

In addition to reaching single-atom sensitivity in the 
readout there is a range of possibilities to extend our 
detection capabilities. One powerful option is to make 
use of the lattice periodicity. A simple averaging over 
lattice sites will reduce the uncertainty in the average 
atom number per site. More detailed information is avail- 
able however, correlating the imaging data of neighbor- 
ing sites. The use of correlation functions has already 
allowed us to separate various noise sources in our mea- 
surements of sub-Poissonian atom number fluctuations. 
Whereas averaging can show that the average number 
of collective excitations within one site is 1, correlation 
functions can be to used show that double occupancy is 
suppressed [3T] , In a similar fashion correlation func- 
tions should boost the sensitivity to detect entanglement 
between neighboring sites. 
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D. Site-selective electric fields IV. SCALING DOWN TO SUB-OPTICAL 

PERIOD LATTICES 



An interesting alternative to optical single-particle de- 
tection is the use of sharp tips for ionization-based detec- 
tion. Arrays of such tips have been fabricated in silicon 
for field emission purposes [12]. Such arrays could be 
integrated on the atom chip and are site-selective by na- 
ture. In combination with the shift register operation 
that we have demonstrated before [4 , a linear array of 
tips would be sufficient to read out a two-dimensional 
lattice. 

An externally applied homogeneous electric field is 
concentrated at the tip of a protrusion; this is the "light- 
ning rod" effect. Therefore, structuring the surface topol- 
ogy of a conductor can be used to tailor the local electric 
field. At a distance r from the center of a tip of a single 
thin, long post (height /i), the external field is enhanced 
by a factor 7 « h/r [32] . 

For the magnetic lattices of several micrometer peri- 
odicity considered in this section, it should be possible 
to achieve useful site-selective electric fields with conven- 
tional semiconductor technology. The arrays of silicon 
field emitters produced and characterized in Ref. [12] 
are very promising in this regard. With external elec- 
tric fields of a few kV/mm, field emission from the tip 
of « 5 /im-high cone-like protrusion was achieved, cor- 
responding to a field enhancement factor f3 > 10 3 . This 
height is well matched to the present magnetic-film struc- 
tures. 

Using silicon technology, a chip can be fabricated with 
surface protrusions at selected sites. To maintain con- 
ceptual simplicity, the protrusions can be located at po- 
sitions where the magnetic film is absent, to avoid large 
changes in the local magnetic-field structure. Application 
of a homogeneous external electric field (generated by 
voltage difference between the chip and a remote counter- 
electrode) will then lead to substantially enhanced elec- 
tric fields near the tip of the protrusions. Because the 
field is static, only modest conductivity of the surface 
is needed to achieve the field enhancement, compatible 
with the present FePt magnetic-film atom chips. 

The highly inhomogeneous electric field will exert a 
force on the atoms (proportional to their polarizability) 
directed towards the tip. For a strong field and narrow 
tip, the atoms will be accelerated and field-ionized when 
the field is sufficiently high. The resulting ions or elec- 
trons are further accelerated by the electric field towards 
the surface and counterelectrode. This can be detected 
at the single-particle level (e.g. by using a channeltron 
or multichannel-plate detector at the counterelectrode), 
and thus appears promising as a non-optical site-selective 
single-particle detector. 



Magnetic trap arrays provide a range of opportuni- 
ties that experiments so far have only started to explore. 
They can provide versatile, tunable model systems for 
the study of strongly interacting quantum many-body 
systems, in potentially new parameter regimes beyond 
what is accessible with optical lattices. In this section 
we will investigate scaling down a two-dimensional trap 
array to sub-optical dimensions, with a lattice parameter 
of 100 nm. 

This would yield a new implementation of the Hub- 
bard model, the prevailing paradigm for quantum gases 
in optical lattices. Scaling down to substantially smaller 
lattice parameters will increase all relevant energy scales, 
opening up new regimes with qualitatively new physical 
phenomena. All key parameters, including the tunneling 
rate and on-site interaction can be dynamically tuned. 

Virtually unlimited design freedom allows for, among 
others, square, triangular, hexagonal (graphene-like) , 
and Kagome lattice geometries. In addition, controlled 
amounts of disorder and designer defects can be built in. 

In this section we consider some of the physics issues 
of this nanoscale scenario. 



A. Scaling and tunability of Hubbard model 
parameters 

The Hubbard model is widely used to describe quan- 
tum gases in optical lattices [33]. It is characterized by 
the tunneling rate J (also called t) and on-site interac- 
tion U. These energy scales can be naturally expressed 
in units of the lattice recoil, Er — (ttH) 2 /2md 2 , where d 
is the lattice constant [34 . A typical value for optical lat- 
tice experiments to date has been d = 425 nm. Lattices 
with period d = 100 nm or even smaller would allow a re- 
alization of the Hubbard model in a novel regime that is 
as yet not accessible by optical lattices. Tunneling rates 
will become very large, as will the on-site interaction. 
This may bring energy scales such as the superexchange 
rate J 2 /U , marking the regime of quantum magnetism, 
within easier reach. 

The result of reducing the lattice parameter from 425 
to 100 nm is an increase of the recoil energy Er from 153 
nK to 2.75 /iK. In Table [T] we compare a few other energy 
scales for the two situations, based on scaling expressions 
for U and J given in Ref. [35] for sinusoidal potentials. As 
a benchmark we calculate U , J, and J 2 /U for the value 
of the lattice depth Vq / Er that corresponds to the super- 
fluid to Mott insulator transition. For a two-dimensional 
square lattice the transition was found experimentally at 
J/U = 0.06 by Spielman et al [36]. For d = 100 nm 
this corresponds to Vo/Er = 6.2. Table [i] shows that all 
energy scales are raised by one to two orders of magni- 
tude, thus relaxing the extreme requirements as posed 
by the optical lattice parameters. The quantity J 2 /U is 
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the energy scale associated with superexchange (quan- 
tum magnetism). It can be seen that scaling down the 
lattice yields a value for J 2 /U of 8 nK, well above the 
lowest reported temperature [37] . 



d 425 nm 100 nm 

Vq/Er 10.4 6.2 
E R 153 nK 2.75 /xK 
U 46 nK 2.2 fiK 
J 2.7 nK 135 nK 
J 2 /U 0.16 nK 8.1 nK 



TABLE I: Comparison of energy scales for an optical lattice 
(period d — 425 nm) and a magnetic lattice with d — 100 nm. 
The value for the lattice depth Vq/Er has been chosen such 
that U/4J = 4.2, corresponding to the Mott insulator transi- 
tion for a 2D square lattice as measured in Ref. [36] . 

A lattice period of 100 nm will thus bring us into an 
exciting new regime. The trap-to-surface distance will be 
scaled down in proportion, because the magnetic field of 
any Fourier component in the surface magnetization with 
wave vector k decays exponentially as B ~ exp(— kz). 
Answering the question of how far we can ultimately scale 
down will therefore depend crucially on our understand- 
ing of the interaction with the surface. The value of 100 
nm appears challenging but feasible. We briefly discuss 
the orders of magnitude of a few different aspects of the 
surface interaction that will play a role in our case. 



B. Optimized design 

The best design for a magnetic film pattern which 
would result in a microtrap lattice with specific trap- 
ping parameters is not always self-evident, most often 
relying on experience from numerous trials. A fast linear 
programming algorithm has been developed to automat- 
ically generate magnetization patterns which provide op- 
timal atom confinement while respecting desired lattice 
symmetries and trap parameters [23] . 

The goal of the optimisation was to create a lattice of 
the desired geometry with the highest possible trap fre- 
quencies, while having equal barrier heights to all neigh- 
bouring sites, using square or triangular geometries. For 
both geometries, the pattern consists of bands of magne- 
tised material, with edge patterns creating the trapping 
potentials. Around the trap positions, the boundary of 
the pattern can be seen as similar to a Z-wire shape. 
Fig. [4^ a) shows the result for a square lattice. 

The analysis was restricted to single-layer magnetiza- 
tion patterns with binary "step-like" magnetization or 
thickness variations. This is also how the film is most eas- 
ily patterned using the lithographic techniques we have 
at hand. 

The optimization algorithm is an adaptation of those 
used to generate electrode patterns for ion trapping [38 . 
Although there are some critical differences between the 
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FIG. 4: (a) Optimized magnetization pattern for a square 
lattice, (b) Square lattice, biased for symmetric barriers and 
trapping minima at z — 100 nm. Drawn in the plane of the 
trapping minima. Contours are lines of equal B, at 0.5 G 
intervals. Darker shades indicate lower fields, i.e. lower po- 
tential, (c) Same pattern biased to yield ID channels along y, 
again drawn in the plane of the trapping minima with 0.5 G 
contour spacing. 
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two systems on a practical level, mathematically the 
analogies are close. 

The film of magnetised material can be written as 
an integral over dipoles. We describe an infinite peri- 
odic lattice of magnetic film with out-of-plane magneti- 
sation M(p), with p = (x,y) the spatial coordinate, by 
a Fourier series. The magnetic potential is then written 
as: 

(j){r) = ^hM e~ knmZ [C nm cos(/c nm • p) + 

n,m (2) 
H~ $>nm s i n (^nm ' P)] •> 

where Mo is the value of the magnetization, h is the 
height of the film, and k nm = nK\ + mK^ are the recip- 
rocal lattice vectors. 

The lattice field is calculated as B(r) = B ext — V0(r), 
where B ext is a uniform external bias field. The Fourier 
series can in practice be truncated, keeping only terms 
with coefficients above some threshold. In general, 
Fourier components with wavelength short compared to 
the desired atom-to-chip distance have little effect on the 
field at the atoms, and can be ignored. The bias field 
still adds tunability to the final magnetic field and thus 
the potential energy landscape. We previously explored 
this feature in the shift register experiment [4]. 

In Fig. plb) we show the distribution of the mag- 
netic field using a bias field chosen to yield field 
minima at z = 100 nm above the surface with equal 
barriers in the x and y directions. The required bias 
field is #ext = (-0.98,-0.39,0.00) mT, the minimum 
field is Bip — 0.76 mT, and the trapping frequencies 
cj/2tt = (2.1,2.0,0.8) MHz. The trap depth is calculated 
as 0.29 mT and the barriers 1.3 mT. For the Hubbard 
model parameters we find U/J = 72, so that the ground 
state is expected to be deep in Mott insulator regime. 

The bias field can also be tuned to yield low barri- 
ers in the y direction, as shown in Fig. Qc). The field 
minima are again at z = 100 nm. The bias field that 
achieves this situation is £? ext = ( — 1.99, —0.04, 0.00) mT, 
the minimum field is Bip = 1.83 mT, and the trapping 
frequencies uj/2tt = (1.8,1.2,0.6) MHz. The trap depth 
is now 0.16 mT and the barriers 0.023 mT. Remarkably 
the number of field minima has been doubled. With a 
period of 50 nm in the y direction the lattice recoil is 
now four times higher. We now find U/J = 1.4, well 
below the Mott insulator transition. 

The pattern in Fig. Qa) can be produced in FePt us- 
ing conventional microfabrication techniques including e- 
beam lithography followed by reactive ion etching. The 
effect of finite resolution is qualitatively equivalent to 
truncating the Fourier series expansion of the magnetic 
scalar potential. Because the contributions of higher 
Fourier modes (corresponding to small features) on the 
magnetic pseudo-potential decay rapidly with distance 
from the surface (Eq. Q), the shape of the potential 
at the trap position is relatively insensitive to the fine 
details of the magnetization pattern. 



Ultimately, magnetic-film atom chips will be subject to 
similar physical limitations that limit the density of mag- 
netic data storage, as given by the super-paramagnetic 
limit. For small enough feature size thermal fluctuations 
will drive magnetic domains across the energy barrier as- 
sociated with reversing the magnetization. 

C. Van der Waals attraction to the surface 

At the length scale of hundreds of nanometers, atomic- 
solid state interactions become particularly important as 
the trapped atoms approach ever more closely the surface 
of the atom chip. The attractive Van der Waals potential 
shifts the trap minima towards the surface and lowers the 
trap frequency. This has to be compensated for by in- 
creasing the original magnetic trapping potential. How- 
ever the steepness of the trap is limited by requirement 
that it be slower than the Larmor precession frequency, 
as otherwise trap loss through Major ana spin flips can 
occur. 

Near the surface the attractive Van der Waals (VdW) 
potential varies rapidly with distance as —Cs/z 3 . Thus 
the stiffness of the magnetic traps must be large enough 
to overcome the VdW force. We estimate that at a dis- 
tance of 100 nm the trap frequency must be greater 
than approximately 500 kHz. Such high frequencies 
are unprecedented for magnetic traps, but in our case 
they will in fact occur quite naturally because the mag- 
netic field gradients close to the surface are extremely 
high. The atoms are then confined on a very small 
length scale: for 500 kHz the harmonic oscillator length 
yK/2rauj = 10.7 nm, or only twice the s- wave scatter- 
ing length. These comparable length scales can give rise 
to interesting new physics, such as confinement-induced 
resonances [39] . 

We consider atoms trapped very close to a surface, 
z < A « 124 nm, therefore we will neglect retardation 
effects. The potential can be modeled as a harmonic 
trap plus the Van der Waals potential, 

V{z)= l -muj\z-z ) 2 -C 3 /z\ (3) 

The attractive VdW potential shifts the potential mini- 
mum towards the surface and lowers the trap frequency. 
The potential minimum occurs at z t where V'(z t ) = 0. 
This cannot be solved analytically but we can obtain 
a lowest order approximation, linear in C3. We write 
zt = Zt(Cs), so z t (Cz = 0) = zo and calculate the lowest 
order in the power series: z t = Zq + Csz f t (0). Taking the 
derivative of the equation V'(zt) = with respect to C3 
and solving for z' t (Cs) we get 

J(n \ 3z t (C 3 ) , v 

Z ^ ) = - m ^z t (C,f-l2C^ (4) 

so that the shift in trap position is 
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We want the shift to be small compared to zq, which can 
be expressed as a condition on the trap frequency, 



'3C 3 



(6) 



Instead of requiring the shift in trap position to be 
small we can also require that the trap frequency remains 
positive, i.e. V"(z t ) > 0. 



V"(zo) 



12- 



.C, 



(7) 



For this to be positive we require uo > 2^3Cs/mzQ 1 
a factor of 2 stricter than Eq. (J6|. A slightly im- 
proved estimate is obtained by expanding V"(zQ + 5z t ) ~ 
V"(z ) + Sz t V"'(zo). This leads to the even slightly 
tighter requirement 



W > ^crit = V 2(1 + VQ) 



'3C 3 



mzX 



2.63, 



'3C 3 



mzX 



(8) 



In order to calculate numerical values we need the value 
of Cs which is given by [40] 



C 3 



3 e r 



1 



16 £ r 



1 



(9) 



Here e r is the relative dielectric constant of the surface, 
A = 124 nm and T/2tt = 6 MHz are the reduced wave- 
length and natural linewidth of the dominant transition 
from the ground state (for Rb). For metals or high-index 
dielectrics such as silicon, the factor (e r — l)/(e r + 1) 
is approximately one. Taking the value 0.85 for sili- 
con, we get C 3 = 1.3 x 10~ 48 J m 3 . For z = 100 nm 
we get cj cr i t /27T = 688 kHz near a silicon surface. At 
the same time the trap frequency uj must remain small 
compared to the Larmor (spin precession) frequency, i.e. 

^crit < W < U L . 




120 140 

z(nm) 



180 200 



FIG. 5: Potential perpendicular to the surface, expressed in 
magnetic field units, for the situation of Fig. [4^c). For the 
upper curve the Van der Waals has been neglected, for the 
lower curve it has been included. 
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FIG. 6: Zoom out of Fig. [5] showing more clearly the potential 
barrier between the trap position and the trap surface. 



D. Tunneling into the surface 

The small atom-surface distance of 100 nm also raises 
the question about tunneling of trapped atoms into the 
surface. Although the magnetic field provides a poten- 
tial barrier, the attractive Van der Waals potential limits 
the height of this barrier by adding a diverging attrac- 
tive — 1/z 3 core to the trapping potential. This could 
cause trap losses due to tunneling into the surface. The 
potential energy is 

V(r)=g F m F ^ B B(r)-^. (10) 

z 6 

The transmission coefficient for tunneling through this 
barrier can be estimated in the WKB approximation as 

T =^(- 2 LJ> ( * ) - E] ' k )- 



The tunneling rate is exponentially suppressed with bar- 
rier height and the magnetic field rises rapidly close to 
the surface, so that the barrier to the surface quickly be- 
comes much higher than the barrier between neighbors. 

Taking E <C V(z) ~ \ib^ we can define a character- 
istic length scale over which T drops by a factor e as 
t = ft -1 = h/y/EmjiB B (x 1/y/B. The proportionality 
constant is h/y/Sm [Lb ~ 1 nm VmT. For a typical case 
tunneling to the surface therefore appears to be negli- 
gible. For example if we take a 25 nm thick FePt film, 
patterned with a 100 nm period lattice, and trapping po- 
sitions at 80 nm from the surface, the above expression 
yields 1.5 x lO" 163 . 



E. Johnson noise 

A well known cause of trap loss is Johnson noise in 
the proximity of a conducting surface [4TJ [42]. In the 



9 



vicinity of a conducting material, thermal currents cause 
magnetic field noise, which can flip the spins of trapped 
atoms. The effect is sometimes also called near-field 
blackbody radiation [4TH44] . 

Using much thinner films than usual in atom chips ex- 
periments, we will benefit from the fact that the film 
thickness is much less than the skin depth at the rele- 
vant frequency for spin flips (in the MHz range). The 
skin depth is S = y/2/ ' hqujcf, with uj the spin flip transi- 
tion frequency and a the conductivity of the metal; for 
gold: (Jau = 45 x 10 6 S/m. At a 5 G bias field, for exam- 
ple, we have uo = 2tt x 3.5 MHz, and thus S « 40 /mi. For 
lattice periods of 100 nm a magnetic film of only 20 nm 
thickness will be adequate. Since our atom chips are typ- 
ically coated by a 50-nm reflective gold coating, this will 
be the dominant source of Johnson noise induced trap 
loss. We can obtain an estimate by scaling earlier ex- 
perimental results [44 or based on theoretical work [43] . 
The first yields 

r= ( 4 +?j) Coidii + d/t)}- 1 . (12) 

Here Co ~ 88 /is (for Cu), d = 100 nm is the distance, 
and t = 50 nm is the film thickness. This yields T _1 = 
0.023 s. 

The lifetime using expressions in Ref. [43] can also be 
estimated: 

/8\ 2 3xl0 22 fu\3 5 2 d 2 _ 

where a film thickness h = 50 nm and distance from the 
film d = 100 nm were inserted. The two estimates are in 
reasonable agreement with each other. 

These results indicate that we may still benefit from 
using lower conductivity material for the conductive coat- 
ing. Another possibility is to avoid the reflective coat- 
ing at the position of the lattice. Atoms could first be 
trapped at a different region of the chip and transported 
to the lattice either magnetically or optically. As for the 
magnetic film, dielectric magnetic materials would be an 
interesting option to explore in the future. 

F. Sub-optical detection with tailored electric fields 



promising alternative that circumvents the above com- 
plications. 

Furthermore, such tailoring of electric fields should be 
compatible with scaling to the 100-nm or sub- wavelength 
regime described in this section. For instance there is a 
great interest in the use of carbon nanotubes for field 
emission (of electrons from the tip) and field ionization 
(of atoms in vacuum near the tip). Experiments demon- 
strating field ionization of cold Rb atoms near carbon 
nanotubes have recently been reported [47] 08] , showing 
explicitly that the low ionization energy of ground-state 
alkalis allows for field ionization at fields that are sub- 
stantially below what is necessary for field emission of 
electrons from the surface. The use of nanotubes is very 
attractive from a conceptual point of view, and may be 
necessary when scaling down to the dimensions consid- 
ered here. 



V. SUMMARY AND CONCLUSION 

After briefly reviewing our previous result, we have in- 
vestigated the possibilities to develop a novel, scalable 
quantum information platform based on arrays of mag- 
netic microtraps defined on a magnetic film atom chip. 
Two approaches have been presented, characterized by 
the lattice parameter. 

The first option, with a lattice parameter of 5 /am is 
a promising system as a lattice of ensemble qubits. The 
interactions can be switchable through the transient ex- 
citation to Rydberg levels. First measurements on the 
interaction of Rydberg atoms with the surface have re- 
vealed the presence of electric fields due to adsorbed Rb 
atoms. Some options to mitigate the effects of these 
adatom fields have been discussed. 

The second option, with a sub-optical lattice parame- 
ter of 100 nm, would open up the physics of atomic lat- 
tices in new parameter regimes. We estimated the corre- 
sponding energy scales in the Hubbard model, which are 
one to two orders of magnitude higher compared to opti- 
cal lattices. We investigated the conditions to keep Van 
der Waals attraction to the surface under control and 
showed that tunneling into the surface will be essentially 
zero for practical cases. 



The experimental achievement of single-site resolution 
for ultracold atoms in trapped optical lattices is currently 
under active development [45 , 46 . These experiments use 
optics with very high numerical aperture (NA) to reach 
sufficient optical imaging resolution in order to extract 
site-resolved information. A great challenge in these tech- 
niques is the required relative positioning stability (at the 
10-nm scale) of the lattice sites and the imaging optics. 

When atom chips are used t o defin e the sites, our sce- 
nario to use electric fields (Sec. Ill D ) for site-specific de- 
tection via field ionization emerges as an attractive and 
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